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The  broken  rotor  bar  is  an  unexpected  fault  and  a  common  cause  of  induction 
motor  failures  that  threaten  the  structural  integrity  of  electric  machines.  In 
this  paper,  a  new  approach  to  a  broken  rotor  bar  diagnosis,  without  slip 
estimation,  based  on  the  envelope  of  the  stator  instantaneous  complex 
apparent  power  (SICAP)  is  proposed.  The  envelope  is  obtained  from  the 
SICAP  modulation  and  then  transferred  to  a  computer  for  monitoring  the 
characteristic  frequency  and  its  amplitude  using  the  Fast  Fourier  Transform 
(FFT).  For  this  purpose,  the  winding  function  approach  (WFA)  is  used  to 
simulate  the  broken  rotor  bar  occurrence  in  a  squirrel  cage  induction  motor 
(SCIM)  fed  on  direct  torque  control  (DTC).  The  obtained  simulation  results 
confirm  the  interest  and  efficiency  of  the  proposed  technique.  Even  when  the 
induction  motor  is  operating  at  the  no-load  level  condition,  the  proposed 
method  is  also  efficient  to  detect  the  broken  rotor  bar  fault  at  low  slip 
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1.  INTRODUCTION 

Induction  motors  are  widely  used  in  automated  production  systems.  Most  of  these  systems  require 
high  performances  control  to  maintain  a  stable  velocity  independent  of  any  load  disturbances  or  model 
uncertainties.  The  SCIM  has  proven  itself  through  the  success  that  prevails  in  the  field  of  speed  variation  and 
through  an  influx  of  research  work  which  is  exceptionally  intended  [1,  2].  The  particular  reasons  for  this 
acquired  confidence  in  the  SCIM  are  based  on  its  intrinsic  qualities  such  as  its  simplicity  of  construction,  its 
mechanical  robustness  and  low  purchase  and  manufacturing  costs. 

However,  asynchronous  machines  are  often  subjected  during  their  operation  to  several  stresses  of 
different  natures  (excessive  heating,  magnetic  weariness  caused  by  the  electromagnetic  forces  and 
environmental  stresses  that  the  rotor  must  undergo  during  its  usual  use).  The  accumulation  of  these 
constraints  causes  defects  in  the  different  parts  of  the  machine  such  as  stator  short-circuits  [3],  eccentricity 
[4],  broken  rotor  bar  [5],  which  generates  ultimately  stops,  leading  to  production  losses.  To  avoid  these 
unscheduled  shutdowns,  the  researchers  haves  been  working  for  several  years  to  develop 
maintenance  techniques. 

It  is  well  known  that  the  widespread  method  to  detect  broken  rotor  bars  is  motor  current  signature 
analysis  (MCSA),  which  is  based  on  signal  processing  tools  such  as  Fast  Fourier  Transforms  (FFT).  It 
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depends  on  detecting  the  sidebands  (l±2ks)f  (s  is  the  rotor  slip,  f  is  the  fundamental  frequency  and  k=l,  2, 
3 . . .  ),  inducing  in  the  stator  current  spectrum  producing  by  geometric  and  magnetic  unbalances  caused  by 
broken  rotor  bars.  Therefore,  these  sideband  frequencies  appearing  around  the  fundamental  frequency  and  the 
evaluation  of  their  amplitudes  can  be  identify  and  used  as  a  reliable  and  accurate  approach  to  diagnose  rotor 
bar  faults. 

However,  the  MCSA  has  some  drawbacks,  where  the  reliability  and  the  accuracy  are  affected  by  the 
diagnosis  of  rotor  asymmetries  at  very  low  slip.  Therefore,  the  sideband  frequencies  become  quite  close  to 
the  fundamental  frequency  [6-8].  In  addition,  in  the  closed-loop  drives,  the  controller  loop  compensates  the 
fault  effect  i.e.  mask  the  fault  effect  [9,  10].  Therefore,  using  stator  current  analysis  causes  a  difficulty  in 
fault  detection,  which  motivated  researchers  attempted  to  use  other  quantities  for  fault  analysis  [11-16]. 

To  handle  these  problems,  other  methods  based  on  the  amplitude  modulation  of  the  three-phase 
stator  current  induced  by  a  broken  rotor  bar  is  used  in  aid  of  diagnostic.  In  fact,  the  rotor  fault  effect  can  be 
localized  in  the  stator  current  envelope  spectrum  at  frequencies  expressed  by  the  components  2ksf  [17-20]. 

This  paper  addresses  these  difficulties  with  an  innovative  method  based  on  the  amplitude 
modulation  of  the  SICAP  is  the  so-called  envelope  (SICAPE)  to  detect  the  broken  rotor  bars  fault  in  a  DTC- 
fed  induction  motor.  A  second-order  low-pass  filter  is  used  to  eliminate  the  high  frequency  generated  by  the 
DTC  inverter.  The  Fast  Fourier  Transform  (FFT)  is  applied  to  the  filtered  envelope  component  in  order  to 
extract  the  fault  indexes  using  the  normalized  band  amplitude  of  the  spectrum  components.  The  simulation  is 
performed  using  Matlab/Simulink.  The  obtained  results  show  that  the  proposed  method  is  able  to  detect  the 
broken  rotor  bars  fault  at  any  load  level  condition  in  the  SCIM  drives  as  well  as  at  very  low  slip  condition. 


2.  MATERIAL  AND  METHODS 

The  broken  bar  fault  in  the  SCIM  is  modeled  based  on  the  winding  function  approach  WFA.  All  the 
space  harmonics  in  the  machine  are  taken  into  account.  The  WFA  approach  predicts  the  performances  of  the 
differential  equations  model.  This  model  refers  to  the  coupled  magnetic  approach  by  treating  the  current  in 
each  rotor  bar  as  an  independent  variable  [21].  The  induction  motor  model  is  considered  with  the  following 
simplifying  assumptions  [22]. 

-  The  magnetic  circuits  are  unsaturated, 

-  The  inter-bar  current  is  neglected, 

-  The  skin  effects  are  neglected, 

-  The  effect  of  slots  is  neglected, 

-  The  distribution  of  the  magnetomotive  force  in  the  air-gap  is  sinusoidal. 

The  rotor  squirrel-cage  modeling  is  based  on  the  equivalent  diagram  of  (Nr-H)  meshes  as  shown  in 
Figure  1.  Each  mesh  is  substituted  by  an  equivalent  circuit  represented  by  a  rotor  bar  and  a  segment  of  end¬ 
ring,  respectively  under  resistive  and  inductive  nature  (Rb,  Re,  Lb,  Le). 

Rb  and  Lb,  epresent  the  rotor  bar  resistance  and  is  inductance.  Re  and  Le  are  the  end-ring  segment 
resistance  and  its  inductance.  The  rotor  voltage  equations  of  the  Nr  loops  can  be  expressed  as  follow 


Where,  (|)r,k  is  the  rotor  flux  crossing  the  rotor  loop  k. 

The  dynamic  mathematical  model  of  the  induction  motor  can  be  written  in  vector  matrix  as 
[V]  =  [R].[1]+^([L][/]) 


(2) 


Where  the  voltage  vector  matrix  [V]  is  given  by 


[V]  =  [m  [V,]] 


[14]  —  [t^sa  ^sb  t?5£.]^ 


[[Vr]  =  [Vrl 
Where  Nr  is  the  number  of  rotor  bars. 


Vre]l 


(IXWr+l) 


=  0 


(3) 
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Figure  1.  The  equivalent  circuit  of  the  rotor  squirrel-cage 


The  current  vector  matrix  [I]  is  composed 

TtH  rr  r  1  rr  m  f[^s]  [^sa  ^sb  ^sc] 


(lxWr+1) 


[i]  =  [[4]  [/.]]- ,  nr 

f  Eri  —  Erl  ‘r2  —  ^rNr  ^eJ(ixWr+l) 

The  global  resistance  matrix[7?]  is  given  as  follow 

rpi  _  wni) 

^  ^  [[0](3,  pni)  [^r] 

Where,  [-/?.,]  is  the  diagonal  matrix  of  the  stator  resistances  phases  of  (m,  m)  dimensions. 
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Rsi,  Rs2  and  Rss  are  the  identical  resistances  of  stator  phases  winding. 
The  resistance  [7JJ  is  a  symmetric  matrix  {Nr+1,  Nr+1)  given  as 
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Where,  R^r  =  2(7?j,  -I-  i?g),  7?j,is  the  rotor  har  resistance  and  R^  is  the  end  ring  resistance. 
The  global  inductance  matrix  [L]  is  as  follow 

_  [^s]  [i^sr](3  .Wr+l)  1 

[[i^rs](  Wr+1.3)  [^r](Wr+l.Wr+l). 


Where  the  inductance  matrix  [L^]  is  given  by 


(4) 


(5) 


(6) 


(7) 


(8) 
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[Ls] 
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Ms 
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0 
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(9) 


Where,  L^i,  1^2  and  are  the  identical  proper  inductances  of  the  stator  coil  and  Mj  is  the  mutual  inductance 
of  the  stator  phases.  The  inductance  matrix  [L^]  is  a  symmetric  matrix  +  1,N^  +  1)  given  by 


[LA  = 


Lrr  (Mj-j-  Lj,)  Myy  ...  (Myy  L  )  ^e{l,Ny+l) 

(A/j-j-  Ltyy  A^J-J-  A7j-j-  Ijq 

I^yy  ^b)  Z/j-j-  (AZj-j-  ^b)  *  * 

:  Myy  (^Myy  “  L  )  *.  (A/j-j-  “  L  )  I 

Tfj)  Z/,. 


-L 


e()Vr+l,l) 


^rVQNr.Nr) 

-L, 


-L, 


S(iV7-  + 1,^7-+ 1) 


(10) 


Where,  L^r  =  ^mr  +  ^(.Li,  +  Lg),  is  the  magnetizing  inductance  of  each  rotor  loop.  Lj,  is  the  rotor  bar 
leakage  inductance,  Lg  is  the  rotor  end  ring  leakage  inductance  and  is  the  mutual  inductance  between 
two  rotor  loops. 

The  mutual  inductance  matrix[Msr],  between  the  stator  windings  and  rotor  loops  of  (3,Ai^-l- 
1)  dimensions,  which  signified  that  the  vector  of  the  rotor  currents  comprise  {k  +  1)  elements 
corresponding  to  the  number  of  rotor  cage  bars  (Ai^),  plus  the  end-ring  loop  (N^  +  1). 


Msirk 

^slrk+1 

■  MsirN^+i 

Ms2r2k 

^s2rk+l 

■  Ms2rNr+l 

Ms3r3k 

^s3rk+l 

■  Mssrpi^+i 

(11) 


The  mechanical  rotor  equations  of  velocity  and  position  are  represented  as  follows 


\J-^  =  Te-fvLir-TL 


,  dt  ^ 


(12) 


Where,  J  is  the  moment  of  inertia,  is  the  load  torque  and  Tg  the  electromagnetic  torque  produced  by  the 
machine  is  obtained  by 


7;=0,5.[/r^M[;] 


(13) 


In  order  to  simulate  the  broken  rotor  bars  failure  using  MATLAB  software,  the  broken  bar  resistance  Ri,  is 
strongly  increased  by  an  additional  resistance  Ri,p  called  the  defect  resistance,  expressed  by  the 
following  matrix. 


-0  0 
0  0 
0  0 
:  0 


0 

0 

LLbPk 

0 


0 

~Lii,p]c 

LLbPk 
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...  0- 
...  0 

0  : 

0  : 


(14) 


The  resistance  matrix  [i?^]  of  the  rotor  squirrel  cage  is  changed  taking  into  account  the  additional  defect 
resistance  matrix the  new  resistance  matrix  is  defined  as  follow 

[Rrd  =  [RA  +  KfI]  (15) 


3.  DIAGNOSIS  OF  ROTOR  FAILURE  IN  DTC  FED  INDUCTION  MOTOR 

DTC  is  an  AC  drive  technology  that  directly  controls  the  torque  and  magnetizing  flux  of  the  motor 
(Figure  2).  This  method  estimates  the  magnetic  flux  and  the  motor  torque  by  measuring  the  supplied  voltages 
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and  currents.  If  the  torque  or  flux  falls  outside  a  predefined  tolerance  range,  the  inverter  transistors  are 
switched  to  the  next  state  so  that  they  return  to  the  range  as  fast  as  possible  [23]. 

3.1.  Stator  flux  behavior 

In  the  {a,P)  reference,  the  stator  flux  can  be  obtained  by  the  following  equation 

Vs  =  -  R,l,)dt  (16) 

Where,  <Ps  is  the  stator  flux  vector,  (p^o  is  the  initial  stator  flux  vector  of  <Ps  and,  is  the  stator  resistance.  If 
the  voltage  drop  due  to  the  stator  resistance  neglected,  (p^  can  be  written  as  follow 

Vs  =  Vso  +  (17) 


Figure  2.  DTC  block  diagram  for  induction  motor  drive. 


The  estimated  flux  is  carried  out  by  the  integration  of  the  stator  voltage.  The  selection  of  the  voltage  vector 
depends  on  the  position  of  stator  flux. 


3.2.  Torque  behavior 

The  electromagnetic  torque  component  is  proportional  to  the  vector  product  between  the  stator  and 
rotor  flux  ipjand  (p^.  The  electromagnetic  torque  equation  is  giving  by 

Te=\PhiVs  (19) 

3.3.  Switehiug  strategy  developmeut 

Based  on  the  hysteresis  state  of  the  torque,  the  flux  and  the  switching  sector  of  the  stator  flux,  which 
is  designated  by  a  (20),  the  DTC  algorithm  (Figure  2)  selects  the  suitable  inverter  voltage  vector  to  apply  to 
the  induction  motor  from  the  Table  1.  The  outputs  of  the  switching  vectors  in  the  different  stator  flux  sectors 
are  parameters  for  the  inverter  switching  devices. 

Where, 

a  =  =  tarr^  (20) 

The  active  switching  vectors  are  ^  Fg,  and  the  zero  switching  vectors  are  Fq,  Fy  as  shown  in  the  Table  1. 
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Table  1.  Switching  sequences  of  inverter  voltage  vectors 


Sector  (S) 

1 
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3 
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no 
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000 
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000 
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000 

111 
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010 

on 

3.4.  Description  of  SICAP  for  both  healthy  and  fanlty  cases 

In  this  section,  a  theoretical  study  from  which  explain  the  effect  of  broken  rotor  bars  in  the  DTC  fed 
SCIM.  In  the  healthy  case,  the  SICAPE  determined  from  the  rms  values  of  current  and  voltage  must  be 
constant  and  contained  only  the  dc  component.  These  rms  values  of  current  and  voltage  provided  by  the  DTC 
fed  induction  motor  considered  balanced  and  sinusoidal  [24].  The  SICAPE  resulting  is  given  by 

s  =  llsll  (21) 

where, 

s(t)  =  VS.Ps.r;  (22) 

/  .271  .27r\ 

=  K.  +  jVp  =  J- Pdc  [Sa  +  +  S,e-^T)  (23) 

and 

/  .271  .271^ 

is  =  ha  +jhp  =  [Sa  +  3  -1-5,6  3  j  (24) 

However,  when  a  broken  rotor  bar  fault  appears,  an  additional  ripple  component  produces  manifesting  on  the 
cyclically  repeating  envelope  (Figure  5)  at  a  rate  equal  to  twice  the  slip  frequency  =  2s f.  The  apparent 
power  s(t)  becomes 


s(t)  =  So  -h  Sfc .  cos(2TT/dfc.  t  ±  (25) 

Where,  Sq  is  the  component  in  the  healthy  case,  and  the  second  term  in  (25)  is  called  s^(t),  represents  the 
additional  ripple  component  besides  the  dc  component,  which  presents  the  broken  rotor  bar  fault  at  the 
disturbance  frequency  [25]. 

p/(t)  =  2]“=i  4  •  cos(2TT/dfe.  t  ±  a^fe)  ^0  for  k  =  l 
(S/CO  =2;“=iSfe.cos(2Ti/dfc.t±asfc)  ^0forki=l 


4.  RESULTS  AND  DISCUSSION 

The  procedures  for  obtaining  the  SICAPE  devoted  to  signatures  analysis  revealing  the  broken  rotor 
bars  can  be  summarized  as  follows; 

In  the  first  step,  a  modulation  of  the  stator  apparent  power  is  necessary  to  obtain  the  so-called 
envelope  (Figure  3(a)  and  3(b)).  This  envelope  resulting  from  the  modulation  of  the  stator  apparent  power  for 
a  DTC  fed  IM  in  the  healthy  and  faulty  cases  with  one,  two  and  three  broken  bars  under  different  load 
conditions  is  shown  in  the  obtained  results  presented  in  Figure  5.  The  second  step,  the  identification  of  the 
envelope  consists  of  extracting  only  the  positive  peak  of  the  stator  apparent  power  (Figure  4(a)  and  4(b)).  The 
third  step,  a  Low-Pass  Filter  for  the  DTC  fed  IM.  The  stator  apparent  power,  which  comprises  the  three- 
phase  stator  currents  (22),  which  comprises  a  high-frequency  component  due  to  the  frequency  of  the  space 
vector  modulation  (SVM)  of  the  DTC  inverter.  As  a  general  rule,  the  fundamental  frequency  of  the  stator 
current  varies  from  0  to  50  Hz  and  the  sampling  frequency  in  our  case  went  up  to  4  kHz,  according  to  the 
hysteresis  bands  of  flux  and  torque.  This  SVM  component  is  removed  from  the  stator  apparent  power  signal 
by  a  second-order  low-pass  filter.  As  a  result,  the  envelope  is  isolated  from  the  apparent  power  without  a 
significant  SVM  component.  Therefore,  the  filtered  envelope  of  the  stator  apparent  power  is  useful  as  an 
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indicator  to  detect  broken  rotor  bars  fault.  A  fast  Fourier  transform  (FFT)  is  applied  on  the  filtered  envelope 
to  detect  an  eventual  broken  rotor  bars  fault  from  which  the  slip  frequencies  appear  at  2ksf  as  shown 
in  Figure  6. 


0  5  10  15  20  25  30 

Time  (s) 


(a) 

Figure  3.  SICAP  and  its  envelope  for  healthy  case 
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and  faulty  case  (b),  under  different  load  conditions 
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Figure  4.  SICAPE  and  its  filtered  spectrum  for  healthy  case  (a)  and  faulty  case 
(b),  under  different  load  conditions 
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Figure  5  shows  the  envelope  and  its  filtered  spectrum  in  both  healthy  and  faulty  cases  divided 
according  to  the  load  level  conditions.  The  simulation  shows  clearly  in  the  healthy  case,  the  filtered  envelope 
is  presented  only  the  continuous  component  under  no  load  and  medium  load  conditions  (Figure  5(a)  and 
5(d)).  On  the  other  hand,  in  the  faulty  case,  the  amplitude  modulation  induced  by  broken  rotor  bars  fault 
appears  clearly  at  no  load  condition  (Figure  5(b)  and  5(c))  and  the  amplitude  modulation  increased  caused  by 
medium  load  condition  (Figure  5(e)  and  5(f)). 

Figure  6  presents  the  spectrum  of  the  filtered  SICAPE.  In  the  healthy  case.  Figure  6(a)  and  6(d) 
show  that  there  is  no  characteristic  frequency  for  both  no  load  and  medium  load  conditions.  In  contrast. 
Figure  6(b)  shows  the  envelope  FFT  is  revealed  clearly  the  characteristic  frequency  /^^with  a  slight 
amplitude  indicating  the  presence  of  a  broken  rotor  bar  fault  in  spite  of  no-load  condition  (very  low  slip).  The 
/dfc  amplitude  is  increased  either  by  the  broken  rotor  bars  number  as  shown  in  Figure  6(b)  and  6(c),  by  the 
load  ratio  as  shown  in  Figure  6(b)  and  6(e),  or  by  the  both  together  as  shown  in  Figure  6(b)  and  6(f). 

The  investigation  on  the  low-frequency  components  at  2ksf  makes  it  possible  that  the  sampling  rate 
of  the  filtered  SICAPE  can  be  decimated  before  FFT  is  performed,  which  greatly  reduces  the  calculation 
time.  Therefore,  this  method  proved  its  effectiveness  and  reliability  in  simplifying  the  task  of  accurate 
detection  of  broken  rotor  bars  under  different  load  conditions  even  at  no-load  condition.  Instead  of  the 
classical  MCSA  method  that  has  demonstrated  its  inability  to  detect  broken  rotor  bars  at  very  low  slip,  and 
fails  in  detecting  faults  of  closed  loop  applications  [6,  10]. 
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Figure  5.  SICAPE  and  its  filtered  spectrum.  No  load  condition;  a)  healthy  case,  h)  IBRB,  c)  3  BRB. 
Medium  load  condition;  d)  healthy  case,  e)  1  BRB,  f)  3  BRB 
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Figure  6.  FFT  of  SICAPE  filtered  spectrum.  No  load  condition;  a)  healthy  case,  h)  IBRB,  c)  3  BRB. 
Medium  load  condition;  d)  healthy  case,  e)  1  BRB,  f)  3  BRB 
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5.  CONCLUSION 

In  the  scope  of  this  work,  a  new  approach  was  proposed  to  diagnosis  the  broken  rotor  bars  fault  in  a 
DTC-fed  induction  motor  using  the  signature  analysis  of  the  SICAP  envelope  (SICAPE).  The  proposed 
method  was  performed  to  detect  the  broken  rotor  bars  fault  for  different  load  levels  in  the  SCIM  drives  as 
well  as  for  no-load  condition.  The  system  model  was  developed  and  implemented  in  Matlab/Simulink 
environment.  The  simulation  was  carried  out  in  order  to  estimate  the  SCIM  performances  under  the  effect  of 
the  broken  rotor  bar.  Based  on  the  obtained  results,  the  following  conclusions  can  be  drawn;  first,  the  closed- 
loop  has  no  influence  on  the  fault  inductor,  despite  the  elimination  of  ripples  on  the  output  quantities  ensured 
by  the  effect  of  the  DTC  control.  Second,  the  proposed  method  allows  generating  broken  rotor  bars  fault 
signatures  using  EFT,  avoiding  the  digital  methods,  which  require  an  expensive  technology.  Third,  the 
analyzed  component  shows  only  the  signatures  of  a  broken  rotor  bar,  which  makes  it  possible  to  avoid 
incorrectly  attributed  to  the  others  electrical  or  mechanical  defects  (stator  or  eccentricity  faults).  Fourth,  the 
proposed  method  is  able  to  detect  the  broken  rotor  bar  at  the  no-load  condition  and  the  characteristic 
frequency  is  clearly  visible  with  a  low  amplitude.  Lastly,  a  fault  severity  ratio  is  defined  by  the  amplitude  of 
the  characteristic  frequency  component,  which  considered  as  a  good  indicator  of  the  state  of  the  SCIM 
related  to  the  number  of  broken  rotor  bars  and  load  level  condition. 
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